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Abstract 

In this work, we study the minimum/stopping distance of array low-density parity-check (LDPC) codes. An array 
LDPC code is a quasi-cyclic LDPC code specified by two integers q and m, where q is an odd prime and m < q. 
In the literature, the minimum/stopping distance of these codes (denoted by d(q,m) and h(q,m), resp.) has been 
thoroughly studied for m < 5. Both exact results, for small values of q and m, and general (i.e., independent of q) 
bounds have been established. For m = 6, the best known minimum distance upper bound, derived by Mittelholzer 
(IEEE Int. Symp. Inf. Theory, Jun./Jul. 2002), is d(q, 6) < 32. In this work, we derive an improved upper bound 
of d(q, 6) < 20 and a new upper bound d(q, 7) < 24 by using a new concept of a template support matrix of a 
codeword/stopping set. The bounds are tight with high probability in the sense that we have not been able to find 
codewords of strictly lower weight for several values of q using a minimum distance probabilistic algorithm. Finally, 
we provide new specific minimum/stopping distance results for m < 7 and low-to-moderate values of q < 79. 

I. Introduction 

In this paper, we consider the array low-density parity-check (LDPC) codes, originally introduced by Fan in [1], 
and their minimum/stopping distance. Array LDPC codes are specified by two integers q and m, where q is an odd 
prime and m < q. Furthermore, in this work, C(q,m) will denote the array LDPC code with parameters q and m, 
and d(q, m) (resp. h(q, m)) its minimum (resp. stopping) distance. 

Since the original work by Fan, several authors have consider the structural properties of these codes (see, e.g., 
[2-8]). For high rate and moderate length, these codes perform well under iterative decoding, and they are also 
well-suited for practical implementation due to their regular structure [9, 10]. 

The minimum distance of these codes was first analyzed by Mittelholzer in [2], where general (i.e., independent 
of q) minimum distance upper bounds for m < 6 were provided. Subsequently, Yang and Helleseth [3] investigated 
the minimum distance of these codes in an algebraic way by first proving that the codes are invariant under a doubly 
transitive group of "affine" permutations. Then, they proved the general lower bound d(q, 4) > 10, for q > 7, on 
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the minimum distance. In [4], the general upper bounds d(q, 4) < 10 and d(q,5) < 12 on the minimum distance 
were proved. Furthermore, by combining these bounds with the results in [3], it follows that d(q,A) = 10 and that 
d(q, 5) is either 10 or 12, for q > 7. In summary, 

6, if to = 3, with equality for q > 5 [3] 

10, if to = 4, with equality for q > 7 [3,4] 

d(q,m) < ^ 12, if to = 5, with exact value either 

10 or 12 for q > 7 [3,4] 

32, if to = 6 [2]. 

The case to = 6 has not been treated in the literature before, except for in the initial work of Mittelholzer [2]. 
In this work, we will consider this case in more detail as well as the case to = 7, both from an experimental point 
of view and by deriving an improved upper bound on d(q 7 6) and a new upper bound on d(q, 7). 

This paper is organized as follows. In Section II, some of the basic notation is introduced and the definition 
of array LDPC codes is given. The concept of a template support matrix is also introduced. In Section III, an 
heuristic is presented that will be used to infer a candidate template support matrix. The heuristic analyzes the 
graphical cycle structure of support matrices of codewords/stopping sets for different values of q, with m fixed. In 
Section IV, we use the template support matrix found in Section III to formally prove the improved upper bound 
d{q 1 6) < 20. Furthermore, in Section V, we present a template support matrix for to = 7, found by using the 
heuristic of Section III, which is used to formally prove the new upper bound d(q, 7) < 24. In Section VI, new 
minimum/stopping distance results are presented for fixed values of to < 7 and q < 79. Finally, in Section VII, we 
draw the conclusions. 



II. Preliminaries 

The array LDPC code C(q, m), with parameters q and to, has length q 2 and can be defined by the parity-check 
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Since the number of ones in each row of the matrix in (1) is q and the number of ones in each column is to, the 
array LDPC codes are (to, q) -regular codes. Furthermore, it is not hard to see that the parity-check matrix in (1) 
has rank qm — m + 1, from which it follows that the dimension of C(q, m) is q 2 — qm + m — 1. 

In [3], a new representation for H(g, m) was introduced. In particular, since each column of the parity-check 
matrix H(q, m) has to blocks and each block is a permutation of (1, 0, 0, . . . , 0, 0) T , we can represent each column 
as a vector of integers between and q — 1, where 



and (-) T denotes the transpose of its argument, i.e., the 1-positions are associated with the integers modulo q. 
Furthermore, it follows from (1) and the integer representation in (2) that any column in an array LDPC code 
parity-check matrix is of the form 



where x and y are integers between and q— 1. Thus, a column can be specified by two integers x and y. Also, note 
that since there are q 2 distinct columns in an array LDPC code parity-check matrix, any pair (x,y) G {0,...,q— l} 2 
specifies a valid column. 

In the following, the support matrix of a codeword/stopping set will be the submatrix of H(q, to) corresponding 
to the support set of the codeword/stopping set, i.e., we keep the columns of H(q, to) whose column indices 
coincide with the support set of the codeword/stopping set. Also, we will use the integer representation in (2) for 
the columns of the submatrix. 

Furthermore, a template support matrix with parameters to, q, w, and q is formally defined as an to x w 
matrix with entries that are functions of q and such that it is the support matrix (possibly column-permuted) of 
a codeword/stopping set of weight/size w of C(q,m) for all q > q . The specific matrix which results when a 
template support matrix is evaluated for a specific value of q is called an instance of the template support matrix. 



In this section, we describe an heuristic which can be used to derive upper bounds on the minimum/stopping 
distance of array LDPC codes. For simplicity, we will only consider the codeword case (the stopping set case is 
similar). The heuristic is a three-step procedure: 

1) In the first step, pairs of codewords ci G C(q\,m) and c 2 G C(q2,m), q\ < q2 and to fixed, where ci and 
c 2 have the same graphical cycle structure (a concept to be defined later), are identified. 

2) The second step is to infer a candidate template support matrix (which may or may not exist) such that the 
instances for q = q\ and q = q 2 are the support matrices (possibly column-permuted) of the two codewords 
ci and c 2 , respectively. We emphasize here that the inferred matrix is only a candidate template support 
matrix, since a formal proof is needed to show that all instances for q > q , for some qo, are in fact valid 
(possibly column-permuted) support matrices. 




(2) 



(x, x + y, x + 2y, . . . , x + (to - l)y) T (mod q) 



(3) 



III. Deriving Upper Bounds on d(q, m) 
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3) The third step is a formal proof that the instances of the candidate template support matrix are indeed valid 
(possibly column-permuted) support matrices of codewords for all possible values of q larger than or equal 
to q . 

Note that all instances of a template support matrix may not have their columns in the order implied by the parity- 
check matrix in (1). This is obviously not important, since the order of the columns in a support matrix is not 
relevant (independent of the order, it will represent the same codeword/stopping set). 

A. First Step: Graphical Cycle Structure 

Note that for the array LDPC codes there exists a subgroup of the automorphism group which is doubly transitive 
[3]. This means that (by definition), for each codeword c £ C(q,m), there exists a codeword p(c) (obtained by 
permuting the coordinates of c according to the permutation p from this subgroup) having the coordinates (pi,P2), 
for any < pi < p 2 < q 2 — 1, in its support set. Thus, it is always possible to permute any codeword (using 
permutations from this subgroup) such that the corresponding support matrix contains the columns (0, 0, 0, ... , 0) T 
and (q— 1, 0, 1, . . . , q — 2) T . This is the case since these columns will always be in the parity-check matrix H(q, m) 
for all valid values of q and m. In particular, the column (0, 0, 0, ... , 0) T (resp. (q — 1, 0, 1, . . . , q — 2) T ) is generated 
by x — and y = (resp. x = q — 1 and y = 1) using the representation in (3). 

As argued above, the support matrix can be regarded as an m x w matrix of integers modulo q, where w is 
the weight of the underlying codeword. From this matrix we can make a bipartite graph, denoted by G^ 1 ^ — 
G(V^' j \E^), for each pair of rows (i, j), i < j. The vertex set yW) partitions into two distinct sets which we 
denote by yM and respectively. Now, for each distinct entry in the ith row of the support matrix we associate 
a node in the vertex set yM. Thus, if there are two (or more) identical entries in the ith row of the support matrix, 
then they will correspond to the same vertex in yW. Similarly, for each distinct entry in the jth row of the support 
matrix we associate a node in the vertex set Furthermore, there will be an edge from a vertex «W g yW to 
a vertex £ Vd) if and only if there exists a column in the support matrix in which the entry corresponding 
to uW appears as the ith element and the entry corresponding to «W appears as the jth element. In the following, 
we will refer to the graphs G^'^ as the support matrix graphs. For convenience, we let v& denote the vertex in 
VW representing the entry (or entries) with value a in the ith row of the support matrix of a codeword. Also, due 
to the automorphism, we will assume that the support matrix of a codeword contains the columns (0, 0, 0, ... , 0) T 
and (g-l,0,l,...,g-2) T . 

Let Ci £ C(qi,m) and c 2 £ C(<72, tn), q 2 > qi, be two distinct minimal codewords of the same Hamming weight, 
where a minimal codeword is a codeword that does not have the support set of a nonzero codeword as a proper 
subset of its own support set. From each of the corresponding support matrices we build the support matrix graphs 
G^'ft for each pair of rows < i < j < m — 1, as outlined above. The graphs corresponding to ci and c 2 

are denoted by Gc^ and Gc^\ respectively. Now, ci and c 2 are said to have the same graphical cycle structure 
(by definition) if and only if the graphs Gcl^ and Gc^\ f° r eacn P a ir (*>i)> nave tne same number of cycles of a 
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Fig. 1. Support matrix graph G' 0,1 ' for the support matrix in (4). The cycle with dashed edges is the cycle 

(4^4 1 ^4 ^4 1 8^4 0) .4 1 ^^ 1 ) .4 1 2 ^4 ^4 1 ^4 6 ) ) * ^ i . 

given length containing the edge (v^}_ 1+i ^ mod q yV^}_ 1+ . ^ mod >) and also the same number of cycles of a given 
length containing the edge (vq\vq^). 

The basic idea is to identify pairs of (minimal) codewords Ci e C(qi,m) and c 2 G C(q2,m), q 2 > qi and m 
fixed, with the same graphical cycle structure, since if they do not have the same graphical cycle structure, then it 
is likely (although not impossible when q\ or q 2 is small) that their support matrices can not be instances of the 
same template support matrix. Then, for a pair of (minimal) codewords with the same graphical cycle structure, 
we would like to infer a template support matrix such that the instances for q = q\ and q = q 2 are the support 
matrices (possibly column-permuted) of the codewords ci and c 2 , respectively. 

Example 1. Consider the case q = 47 and m = 6. Using a computer search, we have found a (minimal) codeword 
of weight 20. The corresponding support matrix is 

r 42 46 5 36 46 37 31 11 5 43 6 37 43 42 36 31 11 6 -i 
43 8 39 4 43 39 32 28 20 32 16 28 24 24 20 16 8 4 

44 1 11 42 9 2 6 4 44 11 42 9 5 6 4 1 5 2 (A\ 
45 2 14 45 14 8 8 27 27 21 37 21 37 33 35 35 33 2 \V 
46 3 17 1 19 14 16 1 3 45 16 18 14 17 19 18 46 45 
L 4 20 4 24 20 24 22 26 22 42 26 46 42 46 3 3 43 43 J 

and the support matrix graph G^ ' 1 ^ (corresponding to the first two rows) is shown in Fig. 1. There is one distinct 
cycle in the graph containing the edge (v^q ,i>o^), namely the cycle 

(Jo) (i) (o) (i) (o) (i) (o) (i) (o) (i) (o)\ (5) 
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(indicated with dashed edges in Fig. 1) of length 10. Furthermore, for the support matrix 

p 54 58 5 48 58 49 43 11 5 55 6 49 55 54 48 43 11 6 -i 
55 8 51 4 55 51 38 34 26 38 22 34 30 30 26 22 8 4 
56 1 11 54 9 2 6 4 56 11 54 9 5 6 4 1 5 2 

57 2 14 57 14 8 8 33 33 27 43 27 43 39 41 41 39 2 \ U / 
58 3 17 1 19 14 16 1 3 57 16 18 14 17 19 18 58 57 
L 4 20 4 24 20 24 28 32 28 48 32 52 48 52 56 56 55 55 J 

corresponding to a (minimal) codeword of weight 20 for q = 59 (and m — 6), the corresponding cycle (also of 
length 10) is 

(AO) (i) (o) (i) (o) (i) (0) (i) (0) (i) (0)\ m 

Thus, we get the same cycle lengths. Continuing with the remaining pairs of rows, — (0, 2), (0, 3), ... , (4, 5), 

we get the same cycle lengths for cycles containing the edge (v^ 1+i ( mod q \, ( mod q \) or the edge {vq\vq^) 

for both support matrices. Thus, we would expect that there might exist a template support matrix whose instances 
(possibly column-permuted) for q = 47 and q — 59 are the support matrices in (4) and (6), respectively. 



B. Second Step: Inferring a Candidate Template Support Matrix 

In this subsection, we consider the second step of the procedure, i.e., to infer a candidate template support 
matrix from two minimal codewords with the same graphical cycle structure. This is done by solving simple 2-by-2 
equation systems and congruences. We remark that this procedure will give a candidate template support matrix, 
since we formally need to prove that the resulting matrix is a template support matrix. 

Now, let 

vt = (v {€) v {€> v {j) v {i) ) (8) 

v l IV.o' "ai,ii •■•> y ai i2i -2' °cti, 2i-i' "ai,2i > yo ' 

(i i) 

denote a cycle of length 21, where a^o = oi\^i, in the support matrix graph Gd computed from a given minimal 
codeword ci G C(qi,m). In a similar manner, we denote by 

v, = ( v {€) v {€> v {j) v {i) ) (9) 

2 * "2,0 ' %,1 ' ' ' ' ' %,2I-2 ' %,21-1 ' 02,21 ' V ^ ' 

where a 2 ,o = oi 2t 2i, cycle of length 21 in the support matrix graph computed from a given minimal codeword 

c 2 G C(<Z2, m), where q 2 > Qi- We assume here that Ci and c 2 have the same Hamming weight and also the same 
graphical cycle structure. Now, the purpose is to infer the entries in a matrix 

X Xi ■■■ Xui-l 

XQ+Vtt xi+j/i ••• 2; ra _i+j/ ro _i 

_x + (m-l)y x 1 + (m-l)y 1 ■■■ x w -i + (m-l)y w -i . 

where w is the Hamming weight of ci and c 2 , and such that the instances of the matrix for q = qi and q = q 2 are 
the support matrices (possibly column-permuted) of Ci and c 2 , respectively. 

Algorithm 1 presents such an algorithm, where ip\(va\) (resp. -02 (^S)) denotes the set of column indices of the 
support matrix of Ci (resp. c 2 ) containing the entry a\ (resp. a 2 ) in the ith row. All entries in the resulting matrix 
(after applying Algorithm 1) should be reduced modulo q to get an instance for a specific value of q. The algorithm 
works on two cycles of the same length, one from the support matrix graph of a minimal codeword ci G C{q\,m) 
and the other from the support matrix graph of a minimal codeword c 2 G C(g 2 ,m), where q 2 > q\. The purpose 
is to fill in the entries in a candidate template support matrix, which initially is filled with erasures denoted by *. 
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Algorithm 1 Template Support Matrix Inference 



1: /* Fill in entries in the candidate template support matrix in (10) based on the support matrices of two (minimal) codewords 
Ci 6 C(qi,m) and c 2 € C(q2,m), q 2 > qi, of the same Hamming weight and with the same graphical cycle structure. 
Input: Row indices i and j, and a pair of cycles (vi, V2) (of the same length 21) as defined in (8) and (9). 
Output: A (partial) candidate template support matrix as defined in (10), and a (partial) permutation tt(-). */ 

2: Assign to 1 all integers in {1, . . . , 1} for some integer I > 1. 

3: for r <- to 21 - 1 do 

4: Find an index pair (a, b) (which is also unique) such that a £ ipi (v^ r ) (1 rpi(vai r+1 ) and b £ rp 2 (v^ r ) ^ ^(vi^r+i )> 

where 7 = 1 and S = j if r is even, and 7 = j and 5 = i if r is odd. 
5: Solve the two systems of equations 

^i 1 ' + IVa^ ( mod 9l) = a l,r 
Xa + $y a (mod (Jl) = Ql, r+ 1 

and 

x( b ] + 72/& 2) (mod 92) = Q2,r 
x[ 2) + 5yj; 2) (mod q 2 ) = a 2 , r +i 

6: Find the integers k x and k y in I that give the simplest solutions (for x and y, modulo qiq 2 ) to the two systems of 
congruences 

x = k x ■ x a 1} (mod qi) 
x = k x ■ xf' (mod q 2 ) 



and 



V = k y -y a (mod gi) 
y = kyy ( b 2) (mod g 2 ) 



if \x\ < \x — qiq 2 \ then 



-1 



else 

x (a; - gig 2 ) ■ fc* 1 
end if 

if < \y-qiQ2\ then 

y <-y- ky 1 
else 

2/ <- (y - 3192) • ky 1 

end if 

if x a — * (and y a = *) then 

:r a <— 5, J/a <— y, 7r(6) <s— a, and go to Step 3. 
else if x a / x or y a / y then 

an inconsistency has occurred. Exit, 
end if 
end for 
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Furthermore, the algorithm also updates a permutation ir(-) which gives the index mapping that should be applied 
to the columns of the support matrix of c 2 to get the instance of the candidate template support matrix for q = q 2 . 
The algorithm is run on pairs of cycles (both containing either the edge [v^_ 1+i ^ mod q yV^}_ 1+ . ( mo d <?)) or ^ e 
edge (vq\vq^) as the left-most or first edge in the cycle) until all entries are filled in. In Step 4 of the algorithm, 
an index pair (o, b) is identified, where a (resp. b) is the index of the column in the support matrix of ci (resp. 
c 2 ) containing ol^I (resp. a^l) as the 7th entry and a^+i ( res P- a 2r+i) as ^ e ^ tn entrv - Later in Step 18 of 
the algorithm, these two indices are used to fill the permutation n (n(b) <— a). Actually, the index pairs (a, b) can 
be computed in a preprocessing stage before the algorithm is even run, since they are available by simple cycle 
analysis. 

In Steps 5 and 6 of the algorithm we determine the entries in column a (of the candidate template support matrix) 
based on the two cycles. In Step 5 we first determine the actual values for x and y modulo qi (denoted by and 
y^a\ respectively) for column a of the support matrix of ci, and then the corresponding values modulo q 2 , now in 
column b (denoted by x^ and y^ 2 \ respectively), of the support matrix of c 2 . Then, in Step 6, we find the simplest 
template solutions, i.e., the simplest expressions for both x and y such that they both evaluate modulo q (for q = q\ 
and q = q 2 ) to the correct values as given by the support matrices of the codewords ci and c 2 , respectively, and 
then the entries for x a and y a are filled in the candidate template support matrix as defined in (10) and as indicated 
in Steps 7 to 18. Note that in Steps 8 and 10 neither the inverse nor the product operation are performed and the 
formal string of three characters; x (or x — q\q 2 in Step 10) (with a specific value inserted for x (or x — q\q 2 
evaluated for a specific value of x in Step 10)), •, and k^ 1 (with a specific value inserted for k x ), is assigned to x. 
A similar comment applies to the assignments in Steps 13 and 15. By the simplest solutions we mean the values for 
k x eland k y G I such that, respectively, max(|fc x |, min(|a;|, \x — q\q 2 \)) and maxflfcjj, min(|y|, \y — qiq 2 \)) areas 
small as possible. We remark here that this is only to find a candidate template support matrix with a nice/compact 
representation, which also makes it easier to prove analytically that all instances (possibly column-permuted) are 
indeed valid support matrices of codewords for all values of q larger than or equal to some q, when q is small. 
In any case, for any practical value of q, regardless of the representation, a simple and fast computer search can 
be used to prove that the candidate template support matrix indeed gives a valid support matrix for all values of 
9o < 1 < q, for some qo. For details, see the proofs of Theorems 1 and 2 in Sections IV and V, respectively. Note 
that if a specific pair of values (k x ,k y ) G I 2 gives a valid template support matrix (which of course requires a 
formal proof), then all pairs (k x ,k y ) € I 2 will give a valid template support matrix, possibly with different values 
for q . Interestingly, however, is the fact that the template support matrices will not necessarily be equivalent, in 
the sense that the instances of the different templates, for a given value of q different from q\ and q 2 , will not 
necessarily give the same codeword from C(q,m). In conclusion, a template support matrix is not unique, even if 
it is based on the same pair of codewords. 

Note that Algorithm 1 does in fact identify a one-to-one mapping (through the permutation 7r(-)) between the 
columns of the support matrices of the codewords ci G C{q\,m) and c 2 G C(q 2 ,m) by matching cycles in the 
corresponding support matrix graphs. Then, the template values for x and y are established by matching columns 
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(and solving equations and congruences independently for each column) through this one-to-one mapping. It is 
in fact this particular one-to-one mapping (as opposed to an arbitrary mapping) that make sure that the resulting 
candidate template support matrix have entries that appear an even number of times (the codeword case) or at least 
two times (the stopping set case) in each row. 

In principle, one type of error condition can occur, i.e., we can exit in Step 20. This happens when a previous 
pair of cycles has determined the entries in column a and then the current pair of cycles gives different values. If 
the algorithm exits in Step 20, we need to start from scratch by considering a different pair of minimal codewords 
Ci e C(qi,m) and c 2 G C(q2,m) of the same Hamming weight and with the same graphical cycle structure, and 
revert (back to erasures) all the entries filled in so far in the candidate template support matrix. 

In Step 5 of Algorithm 1 two systems of equations need to be solved. They have the following solutions: 

= a\,r - j(8 - 7)~ 1 (ai, r +i - ai. r ) (mod qi) 

Va ] = ( s - 7) _1 (ai,r+i - ai,r) (mod q x ) 

4 2) = "2,r - 7(5 - iy 1 {a 2 , r +i - a 2 . r ) (mod q 2 ) 

y { b ] = (S - iy 1 {a2.r+i - a 2 , r ) (mod q 2 ) 

which also gives the rationale behind the assignment to the set of integers 1 in Step 2 of the algorithm, since the 
solutions involve a multiplication by (6 - 7) _1 . 

In Step 6 of Algorithm 1 two systems of congruences need to be solved. They have the following solutions: 

x = k x (xW+q 1 -K(x ( ? ) -xP)) (mod qi q 2 ) (11) 

v = k y (y { a ] + ■ <v { b 2) - v { a ] )) (mod^) (12) 

modulo qiq 2 , where n can be found using the extended Euclidean algorithm which yields integers n and rj such 
that k ■ qi + f] ■ q 2 = gcd(qi,q 2 ) = 1. 

We will illustrate the procedure in Example 2 below. 

Example 2. Consider the two cycles in (5) and (7) for q = 47 and q = 59, respectively. Here, i = and j = 1, 
and K = —5 and rj = 4 (since —5 • 47 + 4 • 59 = 1). For r = (see Step 3 in Algorithm 1), a\. r — ai.o = 46, 
a\, r +i = = 0, «2,r = OL2.0 — 58, Qf2.r+i = OL2.1 = 0, 7 = i = 0, and 5 = j = 1. Since 46 appears in the first 
row and in the second row of the third column (column index 2) of the support matrix in (4), a — 2. Similarly, 
b = 2, since 58 appears in the first row and in the second row of the third column of the support matrix in (6). 
This completes Step 4 of the algorithm, and we get the solutions 

46-0- (1- 0)^(0 -46) (mod 47) = 46 

(1- 0)^(0 -46) (mod 47) = 1 

58-0- (1- 0)^(0 -58) (mod 59) = 58 

(l-0)" 1 (0-58) (mod 59) = 1 



T (i) _ 



(2) 

y 2 = 
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in Step 5, from which we can calculate the following solutions for x and y in Step 6 (with I = m — 1 = 5), using 
(11) and (12), respectively: 



k x /ky 


l 


2 


3 


4 


5 


X 


2772 


2771 


2770 


2769 


2768 


x - q x q 2 


-1 


-2 


-3 


-4 


-5 


y 


1 


2 


3 


4 


5 


y - <7i<72 


-2772 


-2771 


-2770 


-2769 


-2768 



Thus, we can fill in 7r(2) — 2, x 2 — — 1, and y 2 = 1 (corresponding to the values k x = 1 and k y = 1, which give 
the simplest solutions). 

In a similar manner, for instance for r = 3, we get oi, r = 0:1,3 = 28, Oi, r +i = 01,4 = 5, o 2 , r = o 2 ,3 = 34, 
02, r+i = 0:2,4 = 5, 7 = j = 1, and S = i = 0. For this case we have a = 9 and b = 9 (from Step 4), and the 
solutions 



T (l) 


= 28 - 


l-(0- 


1) _1 (5 


- 28) (mod 47) 




= (0- 


1) _1 (5 


- 28) 


(mod 47) = 23 


T (2) 


= 34- 


l-(0- 


1) _1 (5 


- 34) (mod 59) 


(2) 
2/9 


= (0- 


1) _1 (5 


-34) 


(mod 59) = 29 



in Step 5, from which we can calculate the following solutions for x and y in Step 6 (with I = m — 1 = 5), using 
(11) and (12), respectively: 





1 


2 


3 


4 


5 


X 


5 


10 


15 


20 


25 


x - qiq 2 


-2768 


-2763 


-2758 


-2753 


-2748 


y 


1386 


2772 


1385 


2771 


1384 


y - 9132 


-1387 


-1 


-1388 


-2 


-1389 



Thus, we can fill in 7r(9) = 9, Xg = 5, and j/g = — 2 _1 (corresponding to the values k x = 1 and k y = 2, which 
give the simplest solutions). 

Continuing with the rest of the values for r (see Step 3 in Algorithm 1) a total of 10 (the cycle length) columns 
of the candidate template support matrix can be determined. To determine the rest of the entries in the matrix, 
other cycle pairs must be considered. For instance, by looking at the graphs G^°' 2 \ we find the cycles 

( v (0) (2) (0) (2) (0) (2) (0)\ 
y v 46 i v l i v 31 ' v i v ' v 9 i v 46 J 

and 

(JO) (2) (0) (2) (0) (2) (0)\ 

for q = 47 a«ii g = 59, respectively. Choose r = 1, from which we get Oi, r = 01,1 = 1, Oi, r+ i = 01,2 = 31, 
02,r = o 2 ,i = 1, o 2 , r +i = o 2 ,2 = 43, 7 = j = 2, one/ (5 = 2 = 0. For this case we have a — 17 one/ 6=17 (from 
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r 


-5 


-1 


5 


-11 


-1 


-10 


-16 


11 




5 




-4 




6 




-10 


-4 


-5 


-11 


-16 


11 


6 -i 





— 4 





8 


-8 


4 


-4 


-8 


17-2 -1 


9 


2~ 


1 


-7-2~ 


1 


17-2 - 


1 


-15-2 -1 9-2 -1 


2" 1 


2 _1 


-7-2 -1 


-15-2 - 1 


8 


4 





-3 


1 


11 


-5 


9 


2 





6 




4 




-3 




11 




-5 9 


5 


6 


4 


1 


5 


2 





— 2 


2 


14 


— 2 


14 


8 


8 


7-2 -1 


7 


2~ 


1 


-5-2~ 


1 


27-2~ 


1 


-5-2 -1 27-2 -1 


19-2 -1 


23-2 _1 


23- 2 — 1 


19-2 -1 


2 








— 1 


3 


17 


1 


19 


14 


16 


1 




3 




-2 




16 




18 


14 


17 


19 


18 


— 1 


-2 


- 





4 


20 


4 


24 


20 


24 


-3-2 -1 


5 


2~ 


1 


-3-2~ 


1 


37- 2 ~ 


1 


5-2" 1 45-2 -1 


37-2 -1 


45-2 -1 


53-2 -1 


53-2" 1 


— 4 


-4 - 



(13) 



Step 4), and the solutions 

xty = 1 - 2 • (0 - 2)~ 1 (31 - 1) (mod 47) = 31 

= (0 - 2)~ 1 (31 - 1) (mod 47) = 32 
4t = 1 - 2 • (0 - 2) _1 (43 - 1) (mod 59) = 43 
y[ 2 7 } = (0- 2)^(43- 1) (mod 59) = 38 

in Step 5, from which we can calculate the following solutions for x and y in Step 6 (with I = m — 1 = 5), using 
(11) and (12), respectively: 





1 


2 


3 


4 


5 


X 


2757 


2741 


2725 


2709 


2693 


x - qiq 2 


-16 


-32 


-48 


-64 


-80 


y 


1395 


17 


1412 


34 


1429 


y - 9192 


-1378 


-2756 


-1361 


-2739 


-1344 



Thus, we can fill in 7r(17) = 17, Xn = —16, and yn = 17 • 2 _1 (corresponding to the values k x = 1 and k y = 2, 
which give the simplest solutions). Continuing (by considering more cycle pairs) we can determine the rest of the 
columns, and we end up with the candidate template support matrix shown in (13) at the top of the page, where all 
entries should be reduced modulo q to get an instance for a specific value of q. The remaining detailed calculations 
are omitted for brevity. 



C. Third Step: A Formal Proof 

The third step is showing that the candidate template support matrix is indeed a valid template support matrix 
for some parameter q n , i.e., the instances for q > q (possibly column-permuted) are all valid support matrices of 
codewords from C(q, m). In fact it is sufficient (to prove an upper bound on the minimum distance) to show that 
the instances for q > qo (possibly column-permuted) all contain as submatrices valid support matrices of codewords 
from C(q,m). In the case an instance (possibly column-permuted) contains as a proper submatrix a valid support 
matrix of a codeword, the established upper bound is obviously not tight. In particular, we need to show, for any 
value of q > q , for some q , that 

1) all entries in a row occur an even number of times, 

2) all columns in the matrix are in fact valid columns in an array LDPC code parity-check matrix, and 

3) at least two columns are distinct modulo q and appear an odd number of times. 
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Note that the first item above will always be satisfied if Algorithm 1 indeed produces a complete candidate template 
support matrix, since the entries in the matrix are filled in by solving systems of equations and congruences that 
are based on cycles (see the discussion in Section III-B). Furthermore, the second item is also always satisfied, 
since by construction all columns are of the form in (3), for some x and y, and all possible values for x and y 
will give a valid column (see the discussion following (3)). Thus, only the third item above needs to be explicitly 
checked if in fact the candidate template support matrix was produced by Algorithm 1 . 

D. Applicability 

The heuristic outlined above in Sections III-A through III-C is very general and can be applied for any pair of 
values (q,m). However, the difficult part is finding low-weight/small-size candidate codewords/stopping sets for 
different values of q, which is increasingly difficult when m grows, since the minimum/stopping distance increases 
with m. For this we have used the algorithm in [11, 12], and the minimum distance probabilistic algorithm in [13]. 

In this work, we have applied the heuristic for m = 6 and m = 7, but remark that it will easily provide the 
upper bounds d(q, 4) < 10 and d(q, 5) < 12 which can be found in the literature [4]. 

E. Improved Algorithm 

The basic algorithm from Sections III-A and III-B can be improved in the sense of increasing its probability of 
success, i.e., finding a valid template support matrix. The key observation in this respect is that even though the two 
codewords ci E C(qi,m) and c 2 £ C(<z 2 , fn) do not have the same graphical cycle structure, their support matrices 
(possibly column-permuted) may still be instances of the same template matrix. The reason is that different entries 
in the template matrix may reduce to the same value modulo q for different values of q. This typically happens 
when either q\ or q 2 is small. A simple way to deal with such scenarios is by relaxing the condition that Ci and 
c 2 should have exactly the same graphical cycle structure. In particular, it may be sufficient to require that the 
minimum cycle length of all cycles containing the edge (v q % }_ 1+i ^ mod q yV^}_ 1+ . ^ mod ^) and the minimum cycle 
length of all cycles containing the edge (vq\v^) are the same for both support matrix graphs and Gc^\ 

0<i<j<m — 1, and then run Algorithm 1 on such pairs of cycles (which have the same length). 

IV. Upper Bound on d(q, 6) 

By using the heuristic from Section III we have found the candidate template support matrix in (13) shown at 
the top of the previous page. All entries in the matrix should be reduced modulo q. At this stage we emphasis that 
this is a candidate template support matrix, since we need to formally prove that the matrix is a template support 
matrix. In particular, we have used the procedure from Section III-B to infer the matrix in (13) from the codewords 
of Example 1, which have the same graphical cycle structure. Also, in Example 2, some of the columns in the 
matrix in (13) were explicitly determined. The rest of the columns can be determined in a similar manner. Details 
are omitted for brevity. We can now prove the following theorem. 
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Theorem 1. The minimum distance d(q, 6) is upper-bounded by 20 for q > 7. 

Proof: The proof is based on the candidate template support matrix in (13). As explained in Section III-C, there 
are three conditions that need to be checked. Also, if the candidate template support matrix was indeed produced by 
Algorithm 1, only the third condition needs to be explicitly checked. For completeness and for providing a formal 
proof, we will however check all three conditions. Obviously, computing an upper bound on the minimum distance 
from a template support matrix based on a codeword is easy; the upper bound is just the number of columns in the 
matrix. Thus, establishing that the matrix in (13) is a valid template support matrix, in the sense that all instances 
(possibly column-permuted) for q > 7 contain the support matrix of a codeword as a submatrix, establishes the 
upper bound of 20, since there 20 columns in the matrix. 

It is easy to check that each entry in each row of the matrix appears exactly twice, which means that the result 
is true if for any value of q > 7 

2) all columns in the matrix are in fact valid columns in an array LDPC code parity-check matrix, and 

3) at least two columns are distinct modulo q and appear an odd number of times. 

Since all columns in the matrix in (13) are of the form in (3), it follows that they are all valid columns in an 
array LDPC code parity-check matrix (see the discussion following (3)). In particular, the values for x, y for the 
first 6 columns are 



X 





-5 


-1 5 


-11 


-1 


y 





1 


1 3 


3 


5 



For the third part of the proof, we need to show, for any value of q > 7, that there exist (at least two) columns 
in the parity-check matrix which are not identical modulo q and appear an odd number of times. This is simple 
(and very fast) to check by a computer search for any finite value of q that would be of any practical value. It is 
only for large values of q that the theoretical proof below is needed. 

Note that the maximum absolute value of the entries in the first row of the matrix in (13) is 16. Thus, the only 
possibility for repeated columns, when q > 2 • 16 — 32, is for two neighboring columns (with identical entries 
in the first row) to be the same. However, by looking at the third row in the matrix, this possibility can be ruled 
out by requiring that q is larger than twice the maximum absolute value of the entries in the third row, i.e., by 
requiring q > 2 ■ 11 = 22. In summary, it follows that there are no identical columns in the matrix in (13) if 
q > max(32, 22) = 32. Furthermore, for values of 7 < q < 32, it can be checked numerically that there are no 
repeated columns in (13), and the result follows. ■ 

We remark that for q = 7, the matrix in (13) reduces to 

'02653645453640323546" 
03016436510531440314 

04142220644422564152 (\a\ 
0520 5 011001313611620 ■ 
06331502135204035465 
.00464363262165152233. 

We observe that there are indeed some identical columns when q = 7. However, the bound in Theorem 1 is still 
valid, since these columns can just be removed from (14) and we will end up in the valid (but column-permuted) 
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3-2 



-1 



-9-2" 



3 ■ 2 _ 1 1 -7-2 
3 ■ 2 _ 1 2 -5-2 
3 • 2 _ 1 3 -3-2 
3 • 2 _ 1 4 -2" 
3-2 -1 5 
3-2 -1 6 



2" 1 
3 ■ 2 — 1 



-5-2~ 
-3-2" 



3-2 
52~ 



-7-2 - 
-3-2 - 
2" 1 
5-2~ 
9-2~ 
13-2 ~ 



-3 9-2 _ 
-1 5-2~ 
1 11-2" 
3 3 
5 13-2 _ 
7 7-2~ 



-7-2 - 

-3-2" 
-5-2~ 



-5 
-15-2 - 
-5-2 _ 
-5-2 _ 




3-2 _ 
11-2" 



-3-2 


3-2~ 



-2 -9-2" 
3-2 -1 -3 

- 1 - 3 • 2 - 
-2 _1 

3-2 _ 



5-2" 
9-2" 
2 

7-2" 
3-2" 
! 5-2- 



-15-2 - 

-2 

3 • 2 _ 1 
13-2 _ 2 
5 



-7-2" 
: -7-2" 



21-2" 
7 



(15) 



support matrix 

["0 2635460325 4" 
030665514431 

041206425615 nfi'l 
052510336162 "> iu ^ 
063121240346 
.00443215152 3. 

which corresponds to a codeword of weight 12, but the bound d(7, 6) < 20 is of course not tight in this case. In fact, 
we found by exhaustive search that the codeword corresponding to the matrix in (16) is indeed a minimum- weight 
codeword. 



Finally, we remark that the template support matrix in (13) for q = 7, 11, 13, 17, and 19 do not give instances 
with columns in the order as implied by the parity-check matrix in (1). This can easily be seen from the sequence 
of y-values for the matrix in (13), which should be nondecreasing. Furthermore, if two y-values are the same, then 
the corresponding sequence of x-values should be nondecreasing. For q > 19, it can easily be proved that the order 
is always according to (1). However, as argued previously, this is not important (independent of the order, a support 
matrix will represent the same codeword/stopping set). 



V. Upper Bound on d(q, 7) 



For the case m — 7 we have found, using the algorithm from [13], the support matrices 

"0 13 7 8 22 6 18 2 21 18 22 2 14 12 13 20 21 7 20 12 14 6 8 " 

13 18 9 8 20 8 4 2 6 13 2 1 2 9 10 20 10 6 8 2 4 

13 2 9 10 1 10 22 14 10 9 13 1 13 13 14 21 22 10 2 21 

13 3 10 11 2 12 1 20 16 16 20 12 1 2 3 10 11 13 17 19 17 19 

13 4 11 12 3 14 3 3 22 4 12 14 15 22 13 3 11 13 13 15 

13 5 12 13 4 16 5 9 5 7 11 11 3 4 11 12 3 16 5 7 9 11 

13 6 13 14 5 18 7 15 11 14 18 22 11 15 16 1 16 6 22 1 5 7 



T 16 10 11 28 9 24 15 17 9 11 2 17 15 16 26 27 10 26 2 27 24 28" 

16 1 11 12 11 26 22 24 18 20 16 2 1 2 12 13 26 13 24 20 18 22 

16 2 12 13 1 13 28 2 27 1 16 16 17 27 28 13 17 13 12 16 

16 3 13 14 2 15 1 7 9 7 9 15 1 2 3 13 14 16 10 6 6 10 

16 4 14 15 3 17 3 14 16 16 18 15 17 18 28 16 3 3 28 4 

16 5 15 16 4 19 5 21 23 25 27 14 3 4 14 15 3 19 25 21 23 27 

.0 16 6 16 17 5 21 7 28 1 5 7 28 14 18 19 1 19 6 18 14 17 21. 

of (minimal) codewords of weight 24 for q = 23 and q = 29, respectively. Note that in the matrix for q = 23 
(the first matrix) the entries 2 and 8 appear four times in the second row, while in the matrix for q = 29 (the 
second matrix) all entries appear twice in the second row. In the first row, however, all entries appear twice for 
both matrices. Thus, there are several cycles containing either the edge {v^lxyV^) or the edge {vq > \vq'^) in the 
support matrix graph G^' 1 ^ (corresponding to the first matrix), while there is only a single such cycle in the support 
matrix graph G^' 1 ' (corresponding to the second matrix), and as a consequence the codewords Ci and c 2 do not 
have the same graphical cycle structure. Note, however, that the minimum cycle lengths are the same, and this is 
also the case for all the other pairs of graphs and Gc?\ 0<i<j<m — 1. Following the discussion in 
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Section III-E, we may apply Algorithm 1, which infers the candidate template support matrix shown in (15) at the 
top of the previous page. Details are omitted for brevity. 
We can now prove the following theorem. 

Theorem 2. The minimum distance d(q, 7) is upper-bounded by 24 for q > 7. 

Proof: The proof is based on the candidate template support matrix shown in (15) at the top of the previous 
page and is almost identical to the proof of Theorem 1 . In particular, it is easy to check that each entry in each 
row of the matrix appears an even number of times and that all columns in the matrix are in fact valid columns in 
an array LDPC code parity-check matrix (all columns are of the form in (3)). 

For the third part of the proof, we need to show, for any value of q > 7, that there exist columns in the parity- 
check matrix which are not identical modulo q and appear an odd number of times. Again, this is simple (and very 
fast) to check by a computer search for any finite value of q that would be of any practical value. It is only for 
large values of q that the theoretical proof below is needed. 

Now, let the largest absolute value of the entries in the ith row of the matrix in (15) which do not involve a 
multiplication by 2 _1 be denoted A,, and let the largest absolute value of the factor in front of 2 _1 of the remaining 
entries in the ith row be denoted by \n. Since a ■ 2~ 1 (mod q), when a is odd (which is always the case in (15)), 
can be written as (q + a)/2, it follows easily that for a row i where all entries are of the form a or a ■ 2 _1 , different 
template entries can never be the same modulo q when q > 2Xi + in. For the first row this bound is 2 • 5 + 11 = 21, 
and for the third row, this bound is 2-2 + 5 = 9. Thus, looking at the first row, the only possibility for repeated 
columns, when q > 21 (the bound for the first row), is for two neighboring columns (with identical entries in the 
first row) to be the same. However, by looking at the third row in the matrix, this possibility can be ruled out by 
requiring that q > 9 (the bound for the third row). In summary, it follows that there are no identical columns in 
the matrix in (15) if q > max(21, 9) = 21. Furthermore, for values of 7 < q < 21, it can be checked numerically 
that there are no repeated columns in (15), and the result follows. ■ 

VI. Numerical Results 

In addition to the analytic results of Theorems 1 and 2, we have performed a computer search to compute the 
exact values for d(q,m) and h(q,m) for small values of q and to. The results are summarized in Table I, where 
the entries that appear in bold are new results. Results from the literature have also been included with an explicit 
reference. 

For in == 6, we have computed the exact value of d(q, 6) for q < 19. For larger values of q, we have run the 
exhaustive algorithm from [11,12] with an upper weight threshold of 16 without finding any codewords. From 
the upper bound from Theorem 1 and the fact that these codes are even-weight codes, we can conclude that the 
minimum distance, for 23 < q < 79, is either 18 or 20. Furthermore, extensive minimum distance calculations 
using the probabilistic algorithm from [13] for several values of q > 23, indicate that the minimum distance is 
indeed 20 for q > 23, from which it follows that the upper bound from Theorem 1 appears to be tight. 
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TABLE I 

Minimum/Stopping Distance Results for Array LDPC Codes for Different Values of q and m 



q 


d( 


9,7) 


d(g,6) 


h(q,5) 


d(g,5) 


h(q,A) 


d(g,4) 




n 

1 


14 


1 O Ml 
IZ [4 J 


9 


TO Ml 
|4J 


Q r^i 
8 [6] 


Q Ml 
8 |4J 


i i 
1 1 


20 


16 [4] 


10 [6] 


10 [4] 


10 [6] 


i n i/ii 
1U [4J 


1 1 


18 or 20 


14 [4] 


12 


12 [4] 


10 


1 A 1/11 
11) [4J 


l 7 
1 / 


< 


24 


16 


12 


1 O Ml 


10 


1 A 1/11 
11) [4] 


1 Q 
1 y 


< 


20 


18 


12 


1 9 IZll 


10 


1 A 1/11 
IU L4J 




< 


22 


J.O 01 


12 


12 


i n 


1 A 1/11 
IU L4J 


7Q 


< 


24 


j.o 01 


12 


12 


10 


1 A 1/11 
11) |4J 


"X 1 
J 1 


< 


24 


1 ft nr On 

± o or ^ u 


12 


12 


10 


1 A 1/11 


37 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 


41 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 


43 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 


47 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 


53 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 


59 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 


61 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 


67 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 


71 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 


73 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 


79 


< 


24 


18 or 20 


12 


12 


10 


10 [4] 



For m = 7, we have only been able to compute the exact value for q = 7 and 1 1 using the exhaustive algorithm 
from [11,12]. For q = 13, we were able to run the exhaustive algorithm from [11,12] with an upper weight 
threshold of 16. In addition, we found a codeword of weight 20 using the probabilistic algorithm from [13], from 
which (and the fact that the array LDPC codes are even-weight codes) we can conclude that the minimum distance 
is either 18 or 20. For larger values of q, 17 < q < 29, the probabilistic algorithm from [13] has provided the 
estimates (in the form of upper bounds) in Table I. Note that even if the results are formally stated as upper bounds, 
the algorithm from [13] indicates that the estimates are indeed likely to give the exact values, which indicates that 
the bound from Theorem 2 is in fact tight (for instance, q = 17 gives a minimum distance of 24 with very high 
probability). For the high values of q (31 < q < 79), Theorem 2 has provided the estimates (in the form of upper 
bounds). 

VII. Conclusion 

In this paper, the minimum/stopping distance of array LDPC codes have been studied. We have presented an 
improved general (i.e., independent of q) upper bound on the minimum distance for the case m — 6, using the 
concept of a template support matrix of a codeword/stopping set, which significantly improves the currently best 
know bound. The bound appears to be tight with high probability in the sense that we have not found codewords of 
strictly lower weight for several values of q using a minimum distance probabilistic algorithm. In addition, we have 
provided the new upper bound d(q, 7) < 24 which also (from extensive numerical computations) appear to be tight. 
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Finally, we have provided several new specific minimum/stopping distance results for m < 7 and low-to-moderate 
values of q < 79. 
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